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Abstract 

> 
m 

O 

I/"") ■ We present a new selection rule for matrix elements of local spin operators 

o , 

in the S = 1/2 "Haldane-Shastry" model. Based on this rule we extend a re- 
cent exact calculation of the ground-state dynamical spin correlation func- 
T ■ tion 5 aft (n,t) = (0|5 a (n, i)5 6 (0, 0)|0) and its Fourier-transform S ab {Q,E) of 

this model to a finite magnetic field. In zero field, only two-spinon excitations 
contribute to the spectral function; in the (positively) partially-spin-polarized 



^ . case, there are two types of elementary excitations: spinons (AS 2 = ±1/2) 

and magnons (AS Z = — 1). The magnons are divided into left- or right-moving 
branches. The only classes of excited states contributing to the spectral func- 
tions are: (I) two spinons, (II) two spinons + one magnon, (Ilia) two spinons 
+ two magnons (moving in opposite directions), (IHb) one magnon. The 
contributions to the various correlations are: S h : (I); S zz : (I)+(II); : 
(I) + (II) + (III). In the zero- field limit there are no magnons, while in the fully- 
polarized case, there are no spinons. We discuss the relation of the spectral 
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functions to correlations of the Calogero-Sutherland model at coupling A = 2. 
75.10.Jm,71.10.+x 
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I. INTRODUCTION 



In 1988 Haldane and Shastry |2|j3[] independently introduced a S = ~ ID spin model on iV 
sites with an exchange interaction that falls off inversely proportional to the distance between 
the spins. In the past few years this model has proven to be solvable to a remarkable extent 
The simple structure of this model even allowed the authors of ref. [[]J to compute the 
zero magnetic field dynamical structure function (DSF) at zero temperature: 

(GS|5£(f)5S(0|GS) = S ab (m -n,t- f), (1) 

(a, b = x,y, z). In this paper we extend these results to a dynamical groundstate correlation 
function in a nonzero magnetic field. Although a closed expression is not available, we are 
able to identify relevant excitations that contribute to these functions. More specifically: 
when we expand the expression ([I]) for the dynamical structure functions in a basis of 
eigenstates of the Hamiltonian {|^)}: 

(GS\S a m (m } t)S b n (n } t')\GS) = 

£<GS| S a m (t) \u)(u\ S b n (t') |GS) e-W-^-Bo) } (2) 

V 

the set of intermediate states \u) that contribute to the sum is finite and contains only 
states that have very small numbers of elementary excitation added to the groundstate, viz. 
two "spinons" and up to two "magnons" . In zero magnetic field the magnon excitations are 
absent (see ) , whereas in very strong magnetic fields (such that all spins are fully polarized 
in the groundstate) only one magnon participates. The set of intermediate states is small as 
a consequence of a new selection rule for matrix elements of the local spin operators between 
eigenstates of this model. We will present it below. 

A more traditional playground for ID spin chains is the Heisenberg model with nearest 
neighbor exchange (NNE) J7|. It shares its low energy properties with the Haldane- Shastry 
model (HSM). However, for the NNE-model the number of excited states contributing to its 
dynamic structure functions is not bounded; consequently these functions are still unknown 
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in both zero and nonzero magnetic field. It seems again that the characterization of the 
HSM as an "ideal spinon gas" || allows one to push forward our understanding of these ID 
spin chains much further. 

For the remaining part of this presentation, we will continue with a brief reiteration of the 
relevant properties of the Haldane-Shastry model (HSM) in Section [TI[ Section |TTT| describes 
the selection rule. We then use this rule to identify the contributions of intermediate states 



to the different DSFs in Sections [TV] to [VT]. Finally in Section |VII| we conclude with a 
comparison to earlier DSF calculations for the NNE-model. 



II. YANGIAN SYMMETRIES 

The Hamiltonian of the HSM in a magnetic field on N sites is given by: 

H = ^T.d- 2 ^- MPij - 1) + h £ s;, (3) 

where d{n) = sin(^) and Pij is the operator that permutes two spins on sites i and j. 
As described in [|2],|5|,|6|] , in zero magnetic field the spectrum of the HSM consists of large 
degenerate multiplets. Responsible for these degeneracies is a symmetry algebra of @ 
identified in || as the Yangian Y(s/(2)) which is generated by the following two vector 
operators: 

N 

jp = St 

i=l 

Ji = ^WijSi x Sj, (4) 

i<j 

where = and the {zA, i = 1 . . . N are iV equally spaced points on the unit circle. 
Both operators commute with the Hamiltonian. Every degenerate multiplet forms a repre- 
sentation of this algebra and is characterized by a so-called (Yangian) highest weight state 
(YHWS). That state is the only member of that multiplet which is annihilated by Jq and 
Ji . The rest of its multiplet is generated when we act upon the heighest weight state with 
the other members of the Yangian algebra. Every Yangian multiplet has an accompanying 
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Drinfeld polynomial P{u) of order < N in u. The roots of this polynomial (which can be 
chosen to be integers or half-integers in the range . . . N) correspond to the elementary 
excitations of this model: spinons. A Drinfeld polynomial can be represented pictorially as 
a sequence r of iV — 1 zeroes or ones, a generalized occupation number configuration. We tag 
on two zeroes in positions and N. As elucidated in [[]] two ones have to be separated by at 
least one 0. The roots of the Drinfeld polynomial are then located between two consecutive 
zeroes. The locations of the ones, {rrii}, are the so-called rapidities. There are M = N ~^ sp 
of them, where N sp is the number of spinons. To clarify this with an example for N = 10 
sites, let us consider the multiplet characterized by the sequence 101000100. It has three 
rapidities mi = 1, m 2 = 3, m 3 = 7 and four spinons causing roots in the Drinfeld polynomial 
at u — 4|, 5~, 8| and 9|. 

Following |4j] , a sequence of m + 1 zeroes (between ones) is to be interpreted as a single 
orbital filled with m spinons in a symmetric state. Since spinons have spin 1/2 this implies 
that the orbital has total spin y (more formally, this sequence gives an "m-string" of roots 
of the Drinfeld polynomial, which, from the representation theory of the Yangian algebra 
||, constitutes a spin- 1 ^ factor in the representation). In the case of the above example 
multiplet that means that the multiplet has total spin content 1 ® 1 = 2 © 1 © and contains 
3x3 = 9 states. With M ones there are therefore M + 1 orbitals — two of which are empty 
in our example. Notice that the maximum possible S z and S tot a state in the multiplet can 
have is N sp /2. That state (the YHWS) has all its spinons polarized up. 

These spinons have a semionic character as is evidenced by the fact that adding two 
spinons (and we always have to add two at a time to avoid getting a sequence with consec- 
utive ones) to a sequence reduces the number of available orbitals M + 1 = ( N ~^ sp + lj by 
1, as opposed to for bosons and 2 for fermions. All the states in a multiplet characterized 
by a sequence {rrii} have the same energy and momentum given by: 

M 2tt 
P = rrii mod N, in units of — , 
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We can also express these quantities in terms of spinon-variables. If we label the M + 1 
orbitals from right to left by spinon momenta —k <k<k = ^r^- spaced by ^ we get: 

P = ^2 k n kv + mod 27r, 

fc 

E = Y, e(k)n ka + 1 £ - fc') + £?(M, iV), (6) 

fccr fcfc',(Tcr' 

where e(fc) = ^(/Cq — ^ 2 ), ^(fc) = f s (fco — \ k\) and is the number of spinons in the orbital 
with momentum k and spin a. E(M,N) only depends on the total number of sites and 
spinons. We also recognize v s as the spinon velocity at the zone boundary, k , in the 
groundstate. 

The action of the Yangian algebra within a multiplet of states is to rotate the spinons 
individually (rather than all of them through the global SU(2) spin operators). For this 
reason Jq~ and annihilate a YHWS since it has all its spinons fully polarized. Therefore 
the YHWS have also been dubbed Fully Polarized Spinon Gas States (FPSG) [|J. In a 
local spin basis {|{ni, . . . ,um})} where the {rii} are the locations of the down spins, the 
wavefunctions of these FPSG states ^(n 1 , . . . ,n M ) = ip(z ni , . . . , z nM ) are polynomials in 
the {z ni } of degree < N. They can be written as x^o ; where \l/o is the h = groundstate 
Jastrow wavefunction, and x is a polynomial known in the mathematical literature as a Jack 
polynomial. Algorithms for their construction exist ||. 

To get a more physical idea of spinon states consider the following wavefunctions in the 
same basis ||: 

*(ni, . . .,n M \o£i, ■ ■ ■,%) = 

2 M M N sp 

n(^-^ f ) n n n ■ ( ? ) 

i<j i=l i=l j=l 

The {a;} are the locations of localized spinons that can point f or j; the {n,} are the positions 
of down spins (other than those of possible localized spinons pointing down). Notice that 
the wavefunction prevents the {n{\ from coinciding with the spinon sites. We call the 
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complement of the set of spinon sites the condensate. It is a singlet under the action of 
total spin. Furthermore for N sp = eq. (0) represents the exact groundstate wavefunction 
for h = 0. The usefulness of these states is limited by the fact that they are not mutually 
orthogonal and worse, overcomplete. However, based on numerical evidence for up to 12 
spinons, it is clear that the space spanned by states (|7|) with a fixed number N sp of localized 
spinons contains only eigenstates of the Hamiltonian belonging to Yangian multiplets with 
N sp or less spinons. 

The subspace of states that have a fixed number of N sp localized spinons, all polarized, 
has the pleasant property that it only contains eigenstates of H with precisely N sp spinons. 
This is clear from the fact that these loc. spinon wavefunctions — although not eigenstates of 
the Hamiltonian — are annihilated by both Jq" and J± (see Appendix A). This is consistent 
with the fact that fully polarized spinon eigenstates are supposed to be of a polynomial form 
with degree < N in the {z ni } [§J, just like wavefunctions (^). In that same article it is also 
shown that these fully polarized localized spinon states are complete as well and span all 
YHWS. 



In a nonzero magnetic field the term: h J2i Sf in the Hamiltonian (^) will give a Zeeman 
splitting within the Yangian multiplets, although its members remain eigenstates of H . As 
a consequence, for increasing magnetic field, the groundstate will contain more and more 
(fully polarized) spinons. In the thermodynamic limit their number is given by: 

N sv / h , 7T 

— £ = 2a = 1 - Jl - — ; h c =-v s , 8 

where a denotes the groundstate magnetization. For h > h c the groundstate is completely 
ferromagnetic. 

The occupation sequence characterizing the Yangian multiplet that contains the ground- 
state, will have the spinons "condense" into the left- and rightmost orbital (in equal num- 
bers), in accordance with the spinon dispersion (^j relation which assumes a minimum at 
±/c . So for a typical magnetic field below the critical value, the groundstate would be the 
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YHWS of a Yangian multiplet that is described by a sequence like 000010101010000. For 
higher fields the 1010101 pattern shrinks as more spinons go into left and right orbitals. 
These wavefunctions happen to be known analytically Pj: 



2 M 



*k • • • >%) = n (** - m^- M+i . (9) 

i<j i=l 

III. STRUCTURE FUNCTIONS AND THE SELECTION RULE 

The dynamic structure function S ab (m — n,t — t') = (GS| S a (m,t)S b (n,t') |GS) with 
S a (m,t) = exp(—itH)S^exp(itH) measures the response of the system to excitations cre- 
ated by flipping/imposing a certain spin on a particular site n in state |GS) at time t' and 
measuring its effect at time t on site m. Obviously, since H conserves S z , only S h , S zz 
and S + ~ are nonzero. At zero temperature the Fourier transform S(Q, E) of its expansion 
(@) in a basis of eigenstates with energy E v and momentum p v looks like: 

S(Q, E)=Y1 M " S ( E ~ ( E » - E o)) S (Q - (P» - Po)) 

M, a = 27r|(z/|5 a (g)|GS)| 2 , (10) 

with S a (Q) = -^=^2^ =1 exp(—inQ)S^, and |GS) the groundstate of the model. I.e. the 
support of S ab (Q } E) in the (Q, E) plane is zero except when (Q, E) corresponds to the 
excitation-energy and -momentum of a state contained in S a (Q) |GS) Jl^ . 



From numerical evidence up to N = 16 sites it has become clear that there are only an 
unexpectedly small number of nonzero matrix elements M®. To resolve parity and other 
accidental degeneracies between the Yangian multiplets we split these degeneracies by ac- 
tually diagonalizing H + XH3, where H3 is the second integral of the motion for this model 
as presented in P,pT|. The eigenvalues of this operator allowed the Yangian occupations 
sequences to be unambiguously identified. States in the multiplets are partially resolved by 
fixing S z and S tot (a unique resolution of states would be obtained by adding another term 
yuJ • Ji to the Hamiltonian. This would correspond to a basis of states within the Yangian 
multiplet obtained through the Algebraic Bethe Ansatz [|T^]). 
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Now let us denote the eigenstates of this model as |T, /x) where T labels a Yangian 
multiplet through an occupation sequence and fi labels the state within the multiplet. Fur- 
thermore define M rr ' (Sf)^i = (T Sf \T /j,' ) . Then the observation made above implies 
that the matrix M rr (Sf) vanishes if the occupation sequences T and V differ "too much" 
in a sense made precise below. This situation is analogous to an ideal gas, where if O is a one 
body operator: (a\0\/3) = if the occupation number configurations of \a) and \[3) differ 
on more than one orbital. E.g. O = p(x) = J2kk' e l( - k ~ k '^ x c\ck> can add or take out a single 
particle from an orbital in an ideal gas, but in an interacting gas it could add unlimited 
numbers of particle/hole pairs. 

However if T and T' do not differ too much, according to the rule there will always be a 
pair of states \i and \J in either multiplet for which the matrix element is nonzero. 

RULE 1 If 7r(T,m,n) is the total number of ones in positions m through n in Yangian 
occupation sequence T then: M rr (Sf) 7^ iff. \7i(T,m,n) — n(T', m, n)\ < 1 for any 
1 <m < n < N — 1. 

The rule is illustrated in Fig. [1|. A general consequence of this rule is that when we choose 
m = 1 and n = N — 1, it follows that the total number of ones in a sequence can't change 
by more than one, i.e. the total number of spinons can only change by +2,0, or -2. It is 
remarkable that according to the rule this also holds on any corresponding subsequences of 
the occupation number sequences. 

The zero magnetic field DSF has been computed in JTJj. The particular structure func- 
tion computed happened to be S~ + (Q,E) (the others are identical because of rotational 
invariance). This function is governed by excitations present in |GS). Since the zero 
field groundstate contains no spinons the rule tells us that we can only expect excitations 
with or 2 spinons. As |GS) has S z = +1 the former is ruled out and in the states in 
the multiplets with 2 spinons, both must be polarized. This was to be expected since we 
can expand £/" | GS ) in a set of localized spinon wavefunctions (0) containing two polarized 
spinons: 
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(ni, . . .,n M ~i\ |GS) = 
N—i 2 1 — (— ) m 

E 77 i _ y -m %,-,^-i|y + m). (11) 
m =i JV 1 z 

So one of the spinons seems to be sitting on the site on which acted and the other is an 
even number of sites removed from it. 

Since we know the spinon dispersion relation, we can demarcate the support of S h (Q, E) 
in the (Q, _E)-plane in Fig. ^. The main steps of the computation of the matrix elements M+, 
i.e. the weight of the DSF at a point (Q,E) on the plot, are the following: since |GS) 
only contains states with two fully polarized spinons, it must be built out of YHWS. These 
wavefunctions are functionally identical to eigenfunctions of the Calogero- Sutherland model 
at coupling A = 2 of particles moving on a ring. Since both wavefunctions are of a polynomial 
form with degree < N the computation of a sum over sites is identical to taking an integral 
over the ring in the continuum model. The action of in the spin chain is translated into 
a particle destruction operator *ff(x,t) (Sf removes a down spin). So the S + ~(Q,E) DSF 
reduces to the Greens function (GS| W(x,t)^(0,0) |GS) in the Calogero-Sutherland model. 
It can be computed in the thermodynamic limit, in which case it can be mapped unto a 
Gaussian hermitian matrix model correlator. The result is: 

S- + (Q, = ± f , 2 K 2 7/f J \ (12) 
8v s \{v 2 s - v()(v 2 s - vi) ) 

with Q = — 2^~(i>i + V2) and E = e(vx) + e^). The DSF matrix element is parametrized 
by V\,V2, which are quickly identified with the velocities of the two spinons in the excited 
state. The (\1'"'" 1 I') Greens function has been obtained recently at finite iV as well ||13|| . 

For h 7^ the three different structure functions S~ + (Q, E), S ZZ (Q,E) and S + ~(Q,E) 
will not be equal, since |GS) is no longer a singlet. In fact it has S z = S tot = = So, 
where N sp = N sp (h) is given by (§). This difference between the three correlation function 



is also expressed in two additional global SU(2) selection rules |TIJ, which rule out certain 
matrix elements based on the total spin S tot and S z of the final state \x inside the Yangian 
multiplet T. 
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In the first place there is the Wigner-Eckart theorem for vector operators such as Sf, 
which tells us that in order for (r,/i| Sf |GS) to be nonzero the total spin S of \T,fi) must 
satisfy Sq — 1<S<So + 1. Secondly, for any state in a multiplet with N sp spinons, we have 
> S > S z — where equality only holds for the YHWS, which has all its spinons f. Now 
put this together with the fact that Sf raises S z by +1, S~ lowers it by 1, and S? leaves 
it the same. Classifying states according to their S z and S tot as types (i) - (vi), following 
Miiller et al. [|Kj we find the following contributions: 



Sf |GS) contains states with S = So + 1 and AiV" sp = +2 (type (Hi) ). 

S- |GS) contains states with S = S + 1 and AN sp = +2 (type (z) ) or S = S and 



AN sp = +2,0 (type (ii 



\o.,b 



S { |GS) contains states with S — Sq + 1 and AN sp = +2 (type (iv) ) or S = So and 
AN sp = +2, (type (v) a ' b ), or S = S - 1 and AN sp = +2, 0, -2 (type (vi) aAc ). 



Since we have an additional quantum number to label states: N sp , we added Latin super- 
scripts a,b,c to the Roman numerals. All 10 contributions are summarized in table (|). 

We will now investigate all three structure functions individually following these selection 
rules. 

IV. S-+(Q,E) 

Type (Hi): AN sp = +2, AS = +1. 

For the occupation sequence of the groundstate in a given magnetic field (e.g. 
000010101010000) let us label the zeroes in the leftmost orbital as the left spinon condensate 
and the ones in the rightmost orbital as the right spinon condensate. From table | we learn 
that action of Sf on the groundstate only produces states with two more spinons, i.e. one 
less 1. This 1 has to come out of the center • • • 10101 • • • region. We can't take more than a 
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single 1 out of the center region — and stow it into the left or right spinon condensate — since 
this would imply a violation of Rule [I] applied to the center region. Taking out a 1 in the 
center region is equivalent to inserting 2 spinons there. A typical nonzero matrix element 
would be ( 0001001001000, fi\ Sf |0001010101000), with the two spinons residing in orbital 
two and three. 

All this means is that we get a simple two spinon spectrum, as in the zero magnetic field 
case. The only difference is that now the momenta of the spinons can only vary from — ko 
to ko, where k = -^(N — N sp ) decreases with increasing magnetic field as N sp = N sp (h) 
according to eq. (||). The support of S~ + (Q, E) is essentially a squeezed version of Fig. ||. 

As for the weight associated with 2-spinon excitations: the calculation for the zero 
magnetic field case ( |l"2"D carries over without problems. The reason for this is twofold. In the 
first place, the nonzero magnetic field groundstate wavefunction (^]) is of the same Jastrow- 
form as the zero field one, with just an extra phase factor ]!« z n t appended. When we take 
the matrix element, the phase factors from ket and bra part cancel each other. Secondly, 
the excited states are again of the YHWS type (S = N sp ) and have to be polynomials. The 
mapping onto a Calogero-Sutherland model matrix element remains therefore legitimate. 

The contribution of just two-spinon YHWS excitations was to be expected since we can 
expand any fully polarized localized spinon wavefunction with N sp spinons acted upon with 

in terms of a set containing N sp + 2 spinons: 

(S^ ai )(ni, . . . ,n M -i) = 

II T 1 _ Zl \ ^( n ii ■ ■ ■ ,n M -i\ai, ■ ■ .,a Nsp ,i,p). (13) 

Here V is a set of M random sites on the circle excluding the spinon sites (see Appendix 



B). Eq. (|TT]) is a special case of this expansion with V equal to the sites that are an even 
number of steps removed from the site on which the local spin operator acts. We could also 
have realized that the number of spinons can't go up by more than 2 when we consider that 
S + and J i annihilate Sf |GS), indicating that the latter must consist of purely YHWS with 
S z = S + 1! 
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V. S ZZ (Q,E) 

For the DSF S ZZ (Q,E) we find similar simple excitations that contribute, although at 
present not all resulting matrix elements can be computed. From the combined selection 
rules we find three types of excitations: (i), (ii) a and (ii) b , in table |I[ They all have in 
common that AN sp = or +2. This isn't surprising: let us consider the state J± (S- |GS )). 
This state is annihilated by and S + , so it must be YHWS (see Appendix C). Therefore, 
since the first action of doesn't change the number of spinons, Sf |GS) must be a mixture 
of states that contain not more than 25*o + 2 spinons, where 2Sq is the number of spinons in 
the groundstate. 

We will now discuss the individual types and where possible compute the values of the 
matrix elements. 

Type (i): AN sp = +2, AS = +1. 

Having identical selection rules, these states sit in the exact same Yangian- and spin 
multiplets as the type (Hi) states, and therefore they also contain two additional spinons. 
However their S z = S tot — 1, so they are no longer YHWS like the type (Hi) states. Never- 
theless they are related by a simple application of S~: 

|r,A7v sp =+2, st °*= s ° +1 ) = i S- T), (14) 

1 sp ' s z =s 1 V 2 (s + 1 ) 1 " v ' 

where |T) denotes the YHWS of the multiplet with occupation sequence Y of type (Hi). 
This allows us to reduce a type (i) matrix element to one that is a multiple of a type (Hi) 
given in eq. (|i~2|) : 



(r.AY„, -V'^I^IGS) 2 



2(50 + 1) 



2(S +1) 



(r| [s + ,s z ] |gs) 2 

(r|S+|GS) 2 , (15) 



where we used that S + annihilates the groundstate. 



13 



Type (ii) a : AN sp = +2, AS = 0. 

This type of state is a member of the same kind of Yangian multiplet as type (Hi) and 
(i) — i.e. with two extra spinons — but it sits in a spin multiplet that does not contain the 
YHWS. We expect the associated matrix element to be proportional to a type (Hi) matrix 
element as well, but we lack the necessary operator that would step us form this state to 
the YHWS. This operator would have to be some member of the Yangian algebra. 

Type (iif: AN sp = 0, AS = 0. 

Since this type of state has AN sp = the number of ones in its occupation sequence 
must be identical to that in the groundstate. As with the previous three types we can delete 
just a single 1 from the center, leaving behind two spinons pointing up. This 1 then must 
go into either the left or the right spinon condensate. Therefore a typical nonzero matrix 
element would be: (00100|l00100l|00000| 5? |00000|l01010l|00000) where the | just helps to 
draw attention to the center region. Rule |1] rules out any additional ones leaving the center 
region. The additional 1 on the left or right can be interpreted as a magnon with S z = —1. 
The limiting case where the magnon "fuses" with two spinons at a boundary between a 
condensate and the center region gives us the groundstate. 

These states must be YHWS since they have = S = S z , like |GS). This fact allows 
us to calculate the corresponding matrix elements. Since both groundstate and excited 
state are YHWS, a mapping onto the Calogero-Sutherland model is valid. In this case we 
need a groundstate density- density correlator (GS| p(x,t)p(x' ,t') |GS) since Sf measures the 
presence or absence of a down-spin (i.e. a particle in the CS-model). This calculation has 



been done by Altshuler et al. fll4| , in the thermodynamic limit, by studying the repulsion of 
energy levels in a random matrix model under a varying perturbation. The energy levels are 
identified with the positions of the particles and the strength of the perturbation corresponds 
to imaginary time. Their original expression depends on three parameters (called A, Ai and 
A2), the latter two of which are compact, and the first one is unbounded. This is precisely 
what we expect from our selection rule: 2 spinons restricted to the center region with 



14 



momenta in the range —ko ■ ■ ■ + &o an d a magnon that can go off all the way to the right or 
left (i.e. ±00 in the thermodynamic limit). In terms of the velocity, v, of the magnon — with 
dispersion relation (H) — and spinon velocities t>i,t>2 their result is as follows: 



\(v,v 1 ,v 2 \p(Q)\GS)\ 2 
(v - v s )(v + v s ) 



(v — vi) 2 (v — V2) 2 



[(v~v l ) + (v-v 2 )] z . (16) 



For completeness we give the relation between the v's and the A's : 

A = v/v s 
\ \ _ (gi+gg) 

AiA2 - ~^r~ 

(1-A?)(1-A^) = (^) 2 . (17) 

Numerical data for S ZZ (Q,E) can be found in Fig. |3] at values of h close to and h c . 
Fig. [| shows the corresponding support of S ZZ (Q,E) in the {Q,E) plane, as predicted from 
the selection rule, in the thermodynamic limit. We notice that for finite-size systems some 
of the features near the lower boundary are absent. 

VI. S+-(Q,E) 

Finally we discuss the S + ~(Q, E) DSF which is governed by the excitations of types (iv) 
- (vi) c that are present in |GS). 

Type (iv): AN sp = +2, AS = +1. 



This type of state is very similar to types (Hi) and (i), as a matter of fact they all 
reside in identical Yangian- and spin multiplets. Therefore type (iv) states differ from the 
groundstate only by two extra spinons in the center region. They are related to their YHWS 
|T) (of type (Hi)) through: 



I S IOI = S n + l 

r,A7V sp =+2, 



So+1 ) = - 1 (s~) 2 |r). (is) 
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Analogous to the calculation for type (i) states, a matrix element of type (iv) can now easily 
be reduced to one involving type (Hi): 



[rAN sp = + 2, s sz : s 3 ^\S!\GS) 

(T\ {S + ) 2 Sr |GS) 



(4S +2)(2S +2) 



(2S + 1)(S + 1) 



(r|5+|GS) 



(19) 



The last matrix element in this equation has already been computed for the S h (Q, E) DSF. 
However the energy of the corresponding excited state here is shifted by 2h in comparison, 
because of the Zeeman term in the Hamiltonian. 

Types (v) a and {vi) a :AN sp = +2, AS = 0,-1. 



States of types (v ) a and(vi) a contain a two-spinon excitation as well — like types (Hi), (i) 
and (iv). However, since they don't reside in the spin multiplet of their YHWS (^f > S), 
the associated matrix elements are unknown. 

Types (v) b : AN sp = 0, AS = 0. 

Type (v) b states are very similar to those of type (ii) b , they only differ in S z by -1. 
Therefore both contain two excited spinons and a single left- or right moving magnon. They 
are related to each other by: 



r,Av s „=o, 



S -l ' 



|r), 



(20) 



and |T) is the type (ii) b YHWS of the multiplet with occupation sequence V. We can now 
trivially relate the matrix elements of type (ii) and (v) b : 



,r,AN sv =o, 



S Z = S -1 



1 5*j I GS ) 



i 

2S 



(T| [S + ;1 Sr] |GS) 



lrl( r l^|GS; 



(21) 



The last matrix element is listed in eq. ([Tq) 
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Type (vi) b : AN sp = 0, AS = -1. 

These states reside in the same Yangian multiplets as the previous type however they 
are not in the spin multiplets of the YHWS: > S = S - 1. Since we lack the operator 
that steps us up to the YHWS, we aren't able to compute the matrix element. 

Type (vif: AN sp = -2, AS = -1. 

This last class of states has AN sp = —2 and therefore the number of ones in its occupation 
sequence goes up by one compared to the groundstate. The selection Rule [I] only allows 
the extra 1 to go into the left or right spinon condensate. As before we are also allowed 
to take a 1 out of the center region and bring it into the left /right spinon condensate. 
Notice however that the rule forbids both of the ones to go into the same condensate: 
one has to be left moving and the other must be right moving. The result is an excited 
state with two magnons and two spinons. A typical nonzero matrix element would be 
(0100|l0010100l|000l|^ |0000|101010101|0000). The | just helps to guide the eye. Also 
present are YHWS of multiplets from the limiting cases where one of the magnons fuses 
with the two spinons in the center; this leaves a multiplet with nothing but one magnon; 
example: (0000|l0101010l|0010 S^ 0000|l0101010l|0000). This single magnon excitation is 
familiar from the strong field regime. 

Since type (vif states are YHWS (^f = S), as is the groundstate, we can repeat 
the calculation of the matrix elements by a mapping onto the Calogero-Sutherland model. 
Because S~ creates a down spin, it corresponds to a particle creation operator in the CS- 
model. The relevant correlation function is therefore (GS| t)W(x', t') |GS). As in the 
(Wty) case for type (Hi) states, a further mapping onto a Gaussian Hermitean matrix 



model allows one to calculate the Fourier transform of this correlation function [|TjJ. The 
result is parametrized by four variables, two of which are compact: Vi,v 2 , and two are 
non-compact:t> , v'\ 

(v,v',v 1 ,v 2 \S-\GS) 2 oc (22) 
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(v 2 - v 2 )(v' 2 - v 2 )\vi - v 2 \ 
(v 2 — v 2 )(v 2 — v 2 )(v 2 — v 2 )(v 2 — vl){v' 2 — v 2 )(v' 2 — v 2 )' 

Energy and momentum in terms of the v's are given by Q = ^-(v + v ' — \(v\ + t^)) and 
E — ^"(f 2 + v' 2 — \{v\ + v 2 ) — v 2 ). It is obvious that the compact parameters are to be 
identified with the spinon velocities and the non-compact ones with the magnon velocities. 

Since we now know all possible excitations contributing to S + ~(Q,E) we can draw its 
support in Fig. [| for low and high values of h. Fig. ^ show numerical data on S + ~(Q,E) 
for those values of h. Table p] summarizes the selection rules and available information on 
matrix elements. 



VII. COMPARISON TO THE BETHE ANSATZ MODEL 



In 1980 Miiller et al. |T(J did a similar calculation of DSFs for the nearest neighbor 
Heisenberg chain. They identified certain types of states called spin wave continuum states 
(SWC) as carrying the dominant contribution to the DSFs. These SWC states can be 
described in a Bethe Ansatz rapidity language by occupation sequences, just like the Yangian 
multiplets in the HSM model. As it turns out, these SWC states correspond to exactly the 
same rapidity sequences that are favored by our selection rules in the HSM! Although the 
dispersion relations for the BA rapidities are different from those in the HSM, the authors 
find the support of the DSFs in the nearest neighbor model to have essentially the same 
shape as we do in the the inverse exchange case. 

However, in the nearest-neighbor Heisenberg chain there are some "anomalous" states, 
characterized by a change of more than ±2 spinons, which contribute to a lesser degree 
to the structure functions and don't lie within the bounds found by the authors. In the 
HSM these contributions are completely absent and once again we find this model to have a 
surprisingly clean structure. So in this sense the Haldane-Shastry model is an ideal spinon 
gas, whereas in the NNE Heisenberg chain the spinons interact. 

Miiller et al. also gave general rules, based on comparing Clebsch- Gordon coefficients, 
determining which matrix elements will survive in the thermodynamic limit. Their conclu- 



sion is that the only surviving ones have S tot = S z . This means only excitations (Hi) for 
S- + (Q, E) , (u) a and (ii) b for S ZZ (Q, E) and (vi) a , (vi) b and (vi) c for S + ~(Q, E) remain rel- 
evant. (Exceptions are single excitations with Q = 0, since these correspond to S + , S z and 
S~ which give macroscopic contributions, as we can see in the figures) . These contributions 
are trivial to compute. 

VIII. CONCLUSION 

We found a remarkably simple selection rule for nonzero matrix elements of local spin 
operators between eigenstates of the Haldane-Shastry model, which is reminiscent of the 
ideal gas single-particle selection rules. One of the consequences of this rule is that the total 
number of spinons can only change by 0, ±2. Within the occupation sequences this holds 
locally as well. 

In the particular case that one of the states in the matrix element is also the groundstate 
in a magnetic field (i.e. fully polarized spinons, condensed into the left- and rightmost 
orbitals in equal amounts), the general selection rule only allows excitations with no more 
than 2 spinons (the rule applied to the center region) and one left- and one right moving 
magnon (the rule applied to the condensates on the left and right) |16| . This implies that 
the structure functions based on the matrix elements involving these states have a finite 
support in a region dictated by convolving the dispersion relations of these particles (Figs.|j 
and ID- These regions have the same shape as those for the nearest neighbor Heisenberg 
chain. The latter model carries some weight outside these regions as well. Therefore the 
HSM model has a much cleaner spinon structure than the BA model. 

Matrix elements that connect a number of the states in these regions to the groundstate 
through the local action of a spin operator have been presented. However, information is 
lacking on those types that involve states, not in the spin multiplet of the YHWS. This 
is particularly bothersome for types (ii) a , (vi) a and (vi) b since these will survive in the 
thermodynamic limit. Their calculation would allow a full reconstruction of the h ^ DFS 
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for the HSM. A more algebraic treatment involving Yangian operators should provide more 
insight. 

This work was supported in part by NSF Grant No. DMR922407 and an IBM Graduate 
Fellowship Award. 

ACKNOWLEDGMENTS 

One of us (JCT) would like to thank Dr. Zachary Ha for useful discussions and making 
available his preprint prior to publication. 

APPENDIX A: 

We show that localized spinon wavefunctions with N sp spinons all pointing j are necessar- 
ily linear combinations of YHWS with exactly N sp spinons. This follows easily from the fact 
that both S + and J+ annihilate these states. If we write ^(n 1: . . . , n M ) = ip(z m , . . . , z nM ) 
where z n = exp(^p) and 

M N sp 

i/j(w 1: . . .,w M ) = - Wj) 2 Y[wi Yl Yl(wi - z a .), (Al) 

i<j i i=l j=l 

then 

(5 + *)(ni,...,n M -i) 

N N 2 ^_, 

= J2^( n ii ■ ■ -,nM-i,j) = • • .,n M -i,e~^) 

3=1 3=1 
= 1>(z ni ,...,z nM ,0)=Q, (A2) 

(where the J2j was recognized as the zero-mode of a Fourier expansion). And with Wjk = 

z 3 +z k . 

Zj-Zk ' 

(J 1 + ^)(m,...,n M -i) ~ 

M N 

"EE^K-'^-i'*)' ( A3 ) 

3=1 i^rij 
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where we used the fact that 2Sf = 1 — 2 J2j=i $i,nj an d that \& is a symmetric function 
which vanishes when two of its arguments coincide. We now use that when convolving a 
polynomial P(z) of degree less than N (such as ip) with Wjk we have: 

N 

5>yP(zi) = NP{ Zj ) - 2z j P'{z j ) - NP(0), (A4) 
i=i 

(P' is the derivative of P). Since ip has a double zero when two of its arguments coincide 
and it vanishes at z = we have \l/ = 0. 

APPENDIX B: 



In this appendix we want to prove eq. (|13|) . Let us first introduce the following identity 
which holds for any set of distinct complex numbers {uo{\: 

M-l M M M 

n & - *) = e n ^-7 - <**). m 

j=l k=l l(^k) Ul UJk i=l 

The RHS is just the Lagrange interpolation formula applied to the function in z on the LHS! 

Say we want to write S^^/{ ai y as a linear combination of localized spinon wavefunctions 
with two more spinons than ^{ ai }, and all spinons pointing f. We fix one of the additional 
two spinons at i, the site on which acts, i.e. : 

(^*{ ai })( n i>---> n M-i) 

= . . . , n M _i, i\aci, a Nsp ) 

= Yl Op*(fti> • • • ,n M -i\ai, . . . ,a Nsp ,i,p). (B2) 

p{¥=i,{ a k}) 

Using eq. (|7]) we can divide out common factors of (z nk — z ni ) etc. , and we are left with: 

Zi]J(z nj - z { ) = a p Y[(z nj - z p ). (B3) 

Jip 



The result follows when we apply the identity (|B 1|) to this equation with z = Zi and a 



z i Hk zl-7 w here the z^ are randomly chosen distinct sites which don't coincide with the 
localized spinons. 
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APPENDIX C: 



We show that S* |GS) can only have or 2 more spinons than |GS) where |GS) is a 
YHWS groundstate in a nonzero magnetic field with 



(ri!, . . . , n M \ GS) = * (^i, • • • , n M ) 



f-M+l 



(CI) 



i<3 



The proof hinges on the fact that in either of those cases (AN sp = or AN sp = +2) acting 
on Sf |GS) twice with will annihilate that state. Potential AN sp = 4, . . . contributions 
should survive as they are at least 2 levels from the top of their Yangian multiplet. Now 



M 



and 



(2S^o) (m,...,n M ) — | 1 — 2 ^ 5j, ni J \&o (m,...,n M ) , 

3=1 



(2J+S'?*o) (n u ...,n M ^) ~ 

M-l N 

~ E E W n k ,l X 
fc=l /^n fe 

M-l 

ni,...,n M _ 1 ,l) 

3=1 

M-l 

2 E uw^ofai, • • • , n M -i,i)- 
k=i 



(C2) 



(C3) 



The first two terms vanish when we apply the convolution theorem with u^, as in Appendix 
A; only the last term survives. Notice that the {n^}; i = 1 . . . , M — 1 cannot be equal to i 
anymore. This state is trivially annihilated by S + as it vanishes at z — 0. Furthermore: 



((j+) 2 ^ ) (m,.,n M - 2 ) = 



M-2 N M-2 

E E *%>,9 E «w*o(wi, • • • > n M-2, g, i) 

p=l q^i,n p k=l 
M-2 N 

+ E E Wn p , q Wq,iVo(n 1 ,...,nM-2,q,i)- 

P=l ?^»,"p 



(C4) 



With the help of the convolution theorem we can set the first term to zero (we can stick 
in the extra term with q = i, which is missing, at no cost since ^ vanishes when q = i). 
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For the second term we use the identity that lies at the heart of the integrability of the 
Haldane-Shastry type models: WijWjk + WjkWki + w^Wij = — 1. 

((j 1 + ) 2 ^)(n 1 ,...,n M - 2 )= (C5) 

M-2 N 

- (! + W i,n p Wn v ,q + W q>i Wi : n p ) \1> (m ,...,n M -2,q,i) ■ 

The first term is zero since \1/ is a polynomial that vanishes at the origin, whereas terms two 
and three can also be put to zero with the help of the convolution theorem. 
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FIGURES 



;r = 100010010000, n\s?\r' = 100010001000,//} ^ o 



;r = 100010010000,// 



Sf 



r" = ooioooooiooo,/'} = o 



FIG. 1. According to the selection rule, there will be states fi in multiplet T that are con- 
nected to others \j! in T' through a local spinon operator. In T and T" there aren't, e.g. since 
K(r,4,8) -7r(r',4,8)| = 2-0 = 2 >1 

-+ 

Support of S (Q,E)fbrh=0 




Q 



FIG. 2. The shaded region show where S ^(Q, E) is nonzero for h = 0. The top boundary cor- 
responds to excitations with two spinons that have identical momentum; on the bottom boundary 
one of the spinons has fixed momentum ±tt. E is given in units of v s /ir. 
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FIG. 3. S ZZ {Q,E) for small h (a = ff = .05) and large h (a = A) on N = 14 sites. £ is in 
units of v s /tt. 




Q Q 

FIG. 4. The regions where S ZZ (Q,E) is nonzero for low and high h. The area shaded dark 
contains the contributions from the excitations with AN sp = +2 (types ((i) and (ii) a ). The 
excitations with AN sp = — type — can also occupy the lightly shaded area. 
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XX 



XX 



Support of S (Q,E) for low h 




2kg 



2n-2nc2n 



Q 



Support of S (Q,E) for high h 




2toj JI 2n-2na 



Q 



FIG. 5. AS XX (Q, E) = S + -(Q, E)+S~ + (Q, E) vs. (Q, E) for low and high h. The contributions 
of AN sp = +2 (types (iv), (v) a and {vi) a ) live in the dark shaded region. These excitations will 
survive in the limit h — > 0. The AN sp = 0, —2 can also occupy the area shaded light. For h — > h c 
only the 1 magnon contributions survive, as can be seen in Fig. |6|. These high field magnon 
excitations are indicated by the thick lines. 



XX 

S (Q,E)atlowh 



XX 

S (Q,E)athighh 




FIG. 6. S XX (Q,E) for low and high h, from numerical simulations on a 10-site chain. E is in 
units of Vs/tt- 
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TABLES 

TABLE I. List of matrix elements contributing to the DSF. The correlation functions in the last 
column refer to those in the Calogero-Sutherland model as they appeared in the previous sections, 
in eqs. (|l^, 15. [H| 19, 21, ^). The entries marked with a (f) do not survive in the thermodynamic 



limit. For the example occupation sequences we act on a groundstate 000|l0101010l|000, with the | 
delimiting the center region. In the column under "Excitation" 5 denotes a spinon and M. denotes 
a magnon. 



DSF 


Type 


gtot 


AN sp 


Typical contributing 
Yangian multiplet 


Excitation 


Matrix element 


S- + (Q,E) 

(S*=S +1) 


{Hi) 


S = S + l 


AN sp = +2 


000|l0010010l|000 


2S 




S ZZ (Q,E) 

(S Z =S ) 


(i)t 


S = S + 1 
S = So 


AN sp = +2 
AN sp = +2 


000|l0010100l|000 
000|l0010100l|000 


2S 
25 





(ii) b S = S AN sp = 100|l0100100l|000 2S + M (pp) 



S+'(Q,E) 




S = 


So + l 


AN sp = 


+2 


000|l0010010l|000 


25 


(S*=S -1) 


( V )0t 


S 


= So 


AN sp = 


+2 


000|01001010l|000 


25 






s = 


So-1 


AN sp = 


+2 


000|01001010l|000 


25 






s 


= S 


AN sp = 


= 


000|l0010010l|010 


25 + M 




(vi) b 


s = 


So -I 


AN sp = 


= 


100|l0010100l|000 


2S + M 




{vif 


s = 


So-1 


AN sp = 


-2 


100|l0010010l|010 


2S + 2M, 



(S +l)(2So+l) 



29 



